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CDMA MULTI-USER DETECTION WITH A REAL SYMBOL 
CONSTELLATION 

CROSS-REFERENCE TO RELATED APPLICATION 

This application claims the benefit of U.S. 
5 Provisional Patent Application No. 60/254,147, filed 
December 8, 2000, which is incorporated herein by 
reference . 

FIELD OF THE INVENTION 

The present invention relates generally to digital 
10 signal processing, and specifically to methods and 
systems for simultaneous reception and processing of 
signals from multiple transmitters sharing a common 
frequency band using code-division multiple access 
(CDMA) . 

15 BACKGROUND OF THE INVENTION 

CDMA is widely used in cellular and other wireless 
systems for multiplexing communications among multiple 
users on a common frequency band. In direct-sequence 
(DS) CDMA, the transmitter multiplies each user's signal 

20 by a distinct code waveform. The code waveforms of the 
different users have the form of "pseudo-noise" (PN) 
sequences, also known as spreading codes, and are chosen 
so as to be as much as possible mutually orthogonal. The 
receiver receives a composite signal made up of the sum 

25 of all of the users' signals, which overlap in time and 
frequency. To recover the individual signal of a given 
user, the receiver correlates the composite signal with 
that user's distinct code. 

Multiple access interference (MAI) limits the 

30 capacity and performance of DS-CDMA systems. MAI results 
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from the lack of complete orthogonality between the 
spreading codes, due to random time offsets between the 
signals of different users, for example. MAI becomes 
increasingly problematic as the number of interfering 
5 users and their power increases, particularly when there 
are power imbalances among the users. These imbalances 
typically result from differences in distance between the 
users and the base station, as well as from fading and 
other factors that affect signal power. The problems 

10 caused by MAI have led a number of groups to develop 
multi-user detection techniques, in which information 
about multiple users and their respective code waveforms 
is used jointly to better detect each individual user. A 
survey of these techniques is presented by Moshavi in 

15 "Multi-User Detection for DS-CDMA Communications," IEEE 
Communications Magazine (October, 1996) , pages 124-136, 
which is incorporated herein by reference. 

To formulate the multi-user detection problem, we 
consider a CDMA receiver that receives signals from K 

20 users. The received complex-envelope baseband signal 
x (t) is a superposition of the users' individual signal 
waveforms, including channel distortion and additive 
noise n(t) : 

K 

25 x(t) = X^Ct) * ^(t) + n(t) (1) 

k = l 

Here denotes convolution, and h k (t) is the impulse 

response of the kth user channel (including delay, 
attenuation, multipath and filtering effects) . In the 
30 derivation that follows, n(t) is assumed to be Gaussian 
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with mean zero and without correlation between successive 
noise samples. 

The kth user's signal waveform is given by: 

5 &(t) = X *>jc(jn)*Wt - mT b )c k (t) (2) 

m = -oo 

In this equation, b k (m) is the jnth symbol transmitted by 
the kth user, wherein b k is selected from a predefined 
constellation (or alphabet) A. For binary phase shift 

10 keying (BPSK) , which is the most common modulation scheme 
in current cellular systems, b k is a bit selected from A = 
{-1,1}. T b is the symbol duration, and p T (t) is a 
rectangular pulse with unity value in t e [0, T) , and zero 
elsewhere. c k (t) is a complex-valued PN spreading 

15 waveform, with a chip period T c = T b /N, wherein N is the 
spreading gain of the DS modulation: 

c k(t) = f ~ \r%\i) + j>f (i)]p rc (t - it c ) ( 3) 

20 and y$\i) are independently drawn from {-1,1} with 

uniform distribution. 

The multi-user detection problem can be stated as 
follows: Let the received signal x (t) be sampled, and 
assume that the impulse responses h k (t) and the spreading 

25 waveforms c^ft; are known. Now estimate the symbols b k (m) 
from the samples x(t) for all k and m. For the sake of 
simplicity in framing the problem, it is assumed that the 
impulse responses of the channels can be expressed as 
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complex scalars h k , that the user signals are mutually 
synchronized, and that the samples are taken at mid-chip, 
at a sampling rate equal to the chip rate (i.e., N 
samples per symbol) . It is possible to relax some of 
5 these assumptions and still achieve multi-user detection, 
as described hereinbelow. 

The symbols transmitted by the K users at symbol 
interval m can be written as a K x 1 real vector: 

10 h(m) := (jbi(m) , . . . , b K (m)f (4) 

The samples taken during the rath symbol interval can be 
arranged as a N x 1 complex vector x(m). Rearranging 
equations (1) and (2) in similar vector form gives the 
15 following expression: 

x(zn) = S(m)b(m) + n(m) (5) 

Here n(m) is a N x 1 complex vector of noise samples. 
20 S (m) is a complex N x K matrix, whose kth column is a 
signature vector s k (m) representing the "symbol signature" 
of the kth user during the rath symbol interval. The 
elements of the signature vectors s k {m) are samples of the 
following waveform: 

25 

s k (t; m) := h k c k {t) , mT b < t < (m + iyr b (6) 

Given the known elements of x(m) and S (m) , in order 
to solve the multi-user detection problem it is necessary 
30 to find the elements of b(m) that best satisfy equation 
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(5) . An optimal, maximum-likelihood solution to this 
problem was framed by Verdu in "Minimum Probability of 
Error for Asynchronous Gaussian Multiple-Access 
Channels," IEEE Transactions on Information Theory 32:1 
5 (January, 1986) , pages 85-96, which is incorporated 
herein by reference. According to this solution, for 
each m, we find the elements of b within the applicable 
constellation A ({-1,1} for BPSK, as noted above) that 

minimize ||x - Sb|j 2 . (For brevity, here and below, we write 
10 simply b, S and x to represent b(ra), S (m) and x (m) . ) 

Although optimal, this detector is usually impractical, 
°3 since it requires an exhaustive search over all possible 
It b e A K . 

As a less costly (though suboptimal) alternative, 
[yl5 Lupas and Verdu proposed a decorrelating detector in 
" J "Linear Multi-User Detectors for Synchronous 

p Code-Division Multiple Access Channels," IEEE 
/"i Transactions on Information Theory 35:1 (January, 1989), 

M= pages 123-136, which is also incorporated herein by 

tl20 reference. The decorrelating detector finds a vector b 
that minimizes fx - Sbf over b e C K , wherein C is the 

complex plane. The "soft decision" vector b is projected 
onto the constellation A K to arrive at the hard decision 
output b . For BPSK modulation, the elements b k are 
25 given simply by the sign of the real part of the 
corresponding elements b k . Assuming K < N and full 
column-rank S, the soft decision solution is given by: 
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wherein S H is the conjugate transpose of S. 

The decorrelating detector has been found to provide 
substantial performance gains over conventional 
5 single-user detection, with significantly lower 
computational complexity than the maximum-likelihood 
detector proposed previously. It performs well even in 
the presence of substantial power imbalances among the 
users. It still entails a substantial computational 
10 burden, however, in inverting the matrix S H S . 

A number of different solutions have been proposed 
O in order to improve multi-user detection performance. 
~£ For example, Duel-Hallen suggests combining the 

decorrelating detector with decision feedback, in 
f jL5 "Decorrelating Decision-Feedback Multiuser Detector for 
i : J Synchronous Code-Division Multiple-Access Channel," IEEE 
s " Transactions on Communications 41:2 (February, 1993), 

W pages 285-290, which is incorporated herein by reference. 
W This detector generates decisions b k iteratively, working 
ri20 down from the stronger user signals to the weaker 
^ signals. The decisions made with respect to the stronger 

users are used in forming decisions for the weaker ones. 
The decision-feedback detector has been found to 
outperform the original decorrelating detector as long as 
25 the bit-error rate (BER) is not too high, and reliable 
estimates of the channel gains h k are available. 

Various other technique have been suggested in the 
patent literature for canceling interference among the 
signals received from multiple CDMA users. For example, 
30 U.S. Patent 5, 644, 592, to Divsalar et al., whose 
disclosure is incorporated herein by reference, describes 
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a method of parallel interference cancellation, which 
estimates and subtracts out all of the interference for 
each user in parallel. U.S. Patents 5,363,403, 

5,553,062, 5,719,852 and 6,014,373, all to Schilling et 
al . , whose disclosures are also incorporated herein by- 
reference, describe methods for solving the (S H S) _1 matrix 
using a fast polynomial expansion. 
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SUMMARY OF THE INVENTION 

Preferred embodiments of the present invention 
provide simpler, more efficient methods and systems for 
detecting CDMA multi-user symbols, particularly BPSK 
5 symbols, using a real statistical model that inherently 
incorporates prior knowledge that the symbols are 
real-valued. As described above, decorrelat ing detectors 
known in the art calculate complex-valued soft decisions 
( b k ) , and then project these values onto the actual 
10 constellation of the user signals to generate the 
corresponding hard decisions {£> k ). In the case of BPSK 
°iS modulation, however, the constellation consists of the 
) y real values ±1, so that b k must be real, even if the user 

SJ waveforms themselves are complex due to the use of 
J7J15 complex PN signatures. 

''■•4 Therefore, in preferred embodiments of the present 

invention, the soft decisions themselves are constrained 
'"-4 a priori to the real axis. This constraint reduces by 
y ; half the number of unknown parameters that must be solved 

1320 for, by eliminating the need to find the imaginary part 
of E k . Application of the constraint thus simplifies the 
decorrelation calculation, improves detection performance 
and doubles the number of users that can be detected 
simultaneously. This constrained decorrelation approach 
25 is applicable not only to BPSK, but also to other 
modulation schemes that use real-valued constellations. 
It can also be used to extract real signals from complex 
signatures that are introduced by factors other than PN 
spreading waveforms, such as the complex response of 
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multiple receiving antennas used for spatial diversity in 
wireless communications. 

In some preferred embodiments of the present 
invention, the novel real statistical model is applied in 
5 a multi-user decision feedback detector. Like decision 
feedback detectors known in the art, the constrained 
decision feedback detector uses decisions made with 
respect to strong user signals to improve the detection 
of weaker user signals. By using the novel real 
10 statistical model, however, the present detector achieves 
better performance than prior art detectors. 

There is therefore provided method for multi-user 
detection, including: 

receiving a complex input signal due to a 
„ L5 superposition of waveforms encoding symbols in a 
real-valued constellation, which are transmitted 
respectively by a plurality of transmitters in a common 
frequency band; 

sampling the complex input signal at sampling 
=^20 intervals over the duration of an observation period to 
provide a sequence of complex samples; 

processing the sequence of complex samples to 
determine soft decision values corresponding to the 
symbols transmitted by the plurality of the transmitters 
25 in the observation period, while constraining the soft 
decision values to be real values; and 

projecting the soft decision values onto the 
constellation to estimate the transmitted symbols. 

In a preferred embodiment, the waveforms include 
30 code-division multiple access (CDMA) waveforms 

transmitted by the plurality of the transmitters, and the 
symbols transmitted by the transmitters are modulated by 
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respective spreading codes to generate the waveforms. 
Typically, the spreading codes include complex-valued 
codes . 

In a further preferred embodiment, the constellation 
5 of the symbols consists of the values +1 and -1, wherein 
projecting the soft decision values includes taking 
respective signs of the soft decision values in order to 
reach a hard decision with respect to the corresponding 
symbols . 

10 The observation period may have a duration 

substantially equal to a single symbol period, during 
which each of the transmitters transmits a single one of 
the symbols, or it may have a duration during which at 
least some of the transmitters transmit more than a 

15 single one of the symbols. 

Preferably, processing the sequence of complex 
samples includes partitioning the samples into real and 
imaginary parts, and processing the real and imaginary 
parts separately to determine the soft decision values. 

20 Typically, the real and imaginary parts of the samples 
are related to the transmitted symbols by an expression 
having a form x = Sb + n, wherein x is a vector including 
real vector elements corresponding separately to the real 
and imaginary parts of the samples, b is a vector 

25 including real vector elements corresponding to the 
values of the symbols, S is a matrix including columns 
corresponding to respective complex signatures of the 
plurality of the transmitter, the columns including real 
entries corresponding separately to the real and 

30 imaginary parts of the signatures, and n is a vector 
including real vector elements corresponding separately 
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to real and imaginary noise components in the samples, 
and wherein processing the sequence of complex samples 
includes inverting the expression to determine the soft 
decision values of the elements of b. 

Preferably, inverting the expression includes 
finding the real values of the elements of b that 

minimize a norm given by ||x - Sb|| 2 , wherein finding the 
real values of the elements of b includes calculating a 
vector b of the soft decision values so that 



Alternatively, inverting the expression includes 
decomposing S to yield an upper-triangular matrix T that 
satisfies an equation z = Tb + v 1 , wherein z and Vi are 
vectors obtained by applying a unitary transformation to 
x and n , respectively, and finding the real values of 
the elements of b iteratively beginning from a final one 
of the elements so as to solve the equation. Preferably, 
processing the sequence of complex samples includes 
ordering the elements of b in an ascending order of power 
of the waveforms transmitted respectively by the 
transmitters, and ordering the entries in S according to 
the order of the elements in b, so that finding the real 
values iteratively includes finding the real values 
beginning from one of the transmitters having a high 
power relative to the other transmitters. 

There is also provided, in accordance with a 
preferred embodiment of the present invention, a 
multi-user receiver, including: 
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input circuitry, coupled to receive a complex input 
signal due to a superposition of waveforms encoding 
symbols in a real-valued constellation, which are 
transmitted respectively by a plurality of transmitters 
5 in a common frequency band, and to sample the complex 
input signal at sampling intervals over the duration of 
an observation period to provide a sequence of complex 
samples; and 

multi-user detection circuitry, coupled to receive 
10 and process the sequence of complex samples so as to 
determine soft decision values corresponding to the 
n symbols transmitted by the plurality of the transmitters 
in the observation period, while constraining the soft 
decision values to be real values, and to project the 
f j_5 soft decision values onto the constellation in order to 
Id estimate the transmitted symbols. 

The present invention will be more fully understood 
|5 from the following detailed description of the preferred 
J1 embodiments thereof, taken together with the drawings in 
M2 0 which: 
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BRIEF DESCRIPTION OP THE DRAWINGS 

Fig. 1 is a schematic, pictorial illustration of a 
multi-user communication system, in accordance with a 
preferred embodiment of the present invention; 
5 Fig. 2 is a block diagram that schematically 

illustrates a transmitter operated by a user in the 
system of Fig. 1; 

Fig. 3 is a block diagram that schematically 
illustrates a receiver with multi-user detection 
10 capability, in accordance with a preferred embodiment of 
the present invention; and 

Fig. 4 is a flow chart that schematically 
illustrates a method for multi-user detection with 
decision feedback, in accordance with a preferred 
15 embodiment of the present invention. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Fig. 1 is a schematic, pictorial illustration of a 
wireless communication system 20 , typically a cellular 
system, in accordance with a preferred embodiment of the 
5 present invention. A base station 22 transmits downlink 
signals to a plurality of users 24, and receives uplink 
signals from the users in return. The signals are 
modulated using a DS-CDMA scheme, with BPSK modulation, 
as described in the Background of the Invention. In 

10 order to separate and process the signals that it 
receives from multiple users 24, base station 22 employs 
a novel detection scheme using a real statistical model, 
as described in detail hereinbelow. 

Fig. 2 is a block diagram that schematically 

15 illustrates elements of a transmitter 26 operated by 
users 24 in system 20. The operation of this transmitter 
is represented mathematically by equations (2) and (3) 
above. Each symbol b(t) to be transmitted is multiplied 
by a complex PN spreading waveform, with different real 

20 (J) and imaginary (Q) parts C x (t) and C Q (t), respectively. 
While b(t) is constrained to the real values ±1, the 
output of transmitter 26 is a complex-valued signal 
b(t)[c x (t) + jC 0 (t)]e :? ' a, ° t , wherein co 0 is the carrier frequency. 

Fig. 3 is a block diagram that schematically 

25 illustrates a receiver 28 used in base station 22 for 
demodulating and decoding the signals from transmitters 
2 6, in accordance with a preferred embodiment of the 
present invention. The signals are received over the air 
by a radio receiver front end 30, and are then digitized 

30 by an analog/digital converter (ADC) 32. The design and 
operation of these elements are well known in the art, as 
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are those of other elements typically used in base 
station receivers, which would normally be included in 
receiver 28 but are omitted here for the sake of brevity. 
ADC 32 provides an input to a multi-user detection block 
5 38 that is substantially of the form given by equation 
(1). A channel estimator 36, as is known in the art, 
estimates the respective impulse responses h k (t) of the 
user channels for use by block 38 . 

Multi-user detection block 38 receives the 

10 superposed signal of all of users 24, modulated by the 
respective PN complex signatures of the users. It 
processes the superposed signal in order to estimate an 
individual symbol stream b k {m) for each user. The novel 
operation of block 38 is described in detail hereinbelow. 

15 Typically, block 38 comprises a digital signal processor 
circuit, comprising a single dedicated chip or a group of 
chips that are configured to carry out the method 
described below. Alternatively or additionally, some or 
all of the functions of block 38 may be carried out in 

20 software on a suitable general-purpose microprocessor. 
Still further alternatively or additionally, certain 
digital processing functions of block 38 that are 
described below may equivalently be accomplished in the 
analog domain using suitable matched filters, as are 

2 5 known in the art. 

The estimated user symbol streams output by 
multi-user detection block 38 are input to an individual 
data processing block 40 for each user. This block 
performs further data decoding and control channel 

30 processing for each user, as is known in the art. The 
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operation of blocks 40 is outside the scope of the 
present invention. 

Multi-user detection block 38 solves equation (5) , 
above, by constraining the soft decision values B k to the 
real axis. This model permits equation (5) to be 
restated as follows in terms of real values only: 

x(m) = SCm)b(m) + n(zn) (8) 



This equation uses the notation S := 




in which 



and S 3 are the separate real and imaginary parts of S, 
respectively. Similar notation is used for the complex 
input signal x and noise n. Equation (8) constitutes a 
real statistical model, which inherently incorporates 
prior knowledge that the symbols are real. Taking 
b e R K , the solution to min||x - Sb|| 2 is the real-valued 
soft decision vector b given by: 

b = (s^^sFx (9) 

Here S T is the transpose of S. The hard decision output 
of block 38 for each user is simply b k = sign(b k ) . 

Because b is constrained to be real, the number of 
real parameters that must be solved for in equation (8) 
is only half the number used in methods known in the art, 
as exemplified by equation (5) . In other words, it is 
evident from this analysis that, at least where BPSK 
modulation is concerned, methods known in the art use an 
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over-parameterized model for their solution. 

Over-parameterization tends to increase estimation error, 
and the inventors have indeed found that in most cases, 
the constrained decorrelation detector implemented in 
5 block 38 achieves a lower bit-error rate (BER) than do 
decorrelation detectors known in the art. The reduction 
in the number of parameters also reduces the 
computational cost of the detector geometrically, so that 
the cost of block 38 is only about one fourth that of a 
10 conventional decorrelating detector for the same number 
of users. On the other hand, for a given spreading gain 
N, equation (8) gives 2N simultaneous equations to use in 
%Q solving for b, which is twice the number given by 

%1 equation (5) . As a result, the usual constraint in 
Sp-5 decorrelating detectors that the number of users K cannot 
y exceed N is now relaxed to K < 2N , permitting receiver 28 

'" 4 to handle twice the number of users as could be detected 

p by a conventional decorrelating detector. 

J1 Although for the sake of simplicity, equations (8) 

M2 0 and (9) relate to samples x(t) collected during a single 
J 15 : symbol interval, the method embodied in these equations 

can readily be extended to multiple successive symbol 
intervals, i.e., to samples collected over [0 , MT b ) . Such 
an extended observation window is typically necessary in 
25 order to deal with loss of synchronism among the user 
signals and to deal with impulse responses h k (t) whose 
duration is longer than a single chip period (due to 
multi-path effects, for example) . In this case, the 
model of equations (1) and (2) can be recast in the 
30 following form: 
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K oo 

x(t) = £ b k (m)s k (t, m) + n(t) (10) 
Jt = ljn = -oo 

Here the symbol signature waveform s k (t;m) of equation (6) 
is replaced by a composite symbol signature waveform, 
5 given by: 

s k (t, m) = h k (t) * [c k (t)p Tb (t - mT b }] (11) 

The composite symbol signature waveform can vary from 
10 symbol to symbol and is of finite duration, say {N+L) 

chip intervals. The size of L reflects the duration of 

the impulse responses. 

If we now restrict our attention to the particular 

observation window [0,MT b ), equation (10) can be rewritten 
15 as follows: 

K 

x(t) = ]T X b k {m)s k {t; m) + n{t) , 0 < t < MT b (12) 
k=l m 

wherein s k (t; m) := s k (t; m)p MTb (t) , and the summation in m is 
20 over all the symbols having non-zero truncated signatures 
s k (t;m) in the observation window. Ordering the symbols 
according to their signature-start-times, the sample 
vector x can now be expressed in the form given by 
equation (5) , x = Sb + n, except that now x is a NM x 1 
25 vector, and the columns of S are the sampled truncated 
signatures. Typically, the columns of S contain a few 
leading zero entries, followed by (N+L) non-zero entries, 
and ending with tailing zeros. 
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Because equation (12) is formally identical to 
equation (5), multi-user detection block 38 can operate 
on asynchronous signals in the same manner as was 
described above with respect to synchronous signals. In 
5 other words, the restriction of b to real values and the 
application of equation (9) to find the values of b k can 
be performed on asynchronous signals, as well. The 
performance of block 38 in the asynchronous case depends 
on the width of the observation window. For ideal 
10 channels and an observation window synchronized with a 
symbol interval of one of the users, b should typically 
O have (2K-1) elements, compared with K elements in the 

"°!f ; synchronous case. For improved performance under 

M= non-ideal conditions, the observation window may be even 

? il5 wider. 

W Fig. 4 is a flow chart that schematically 

illustrates a method for multi-user detection based on 
O decision feedback with the constraint that b is real, in 
ly accordance with a preferred embodiment of the present 
^20 invention. This method can be implemented by multi-user 
detection block 38 in place of the decorrelat ion method 
described above, and gives superior performance when the 
BER is sufficiently low and reliable estimates of the 
channel gains are available. 
25 The method begins by recasting equation (8) in terms 

of an upper triangular matrix T, using QR factorization 
of matrix S , as is known in the art, at a signal 
recasting step 50. The factorization is preferably 
performed using successive Givens rotations or 
30 Householder transformations of the matrix, as described, 
for example, by Golub and Van Loan in Matrix Computations 
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(Johns Hopkins Series in Mathematical Sciences, 1996) , 
which is incorporated herein by reference. Prior to 
factorization, the columns of the matrix are preferably 
ordered according to the approximate strengths of the 
signals received from the respective users, i.e., so that 
the first column corresponds to the weakest user, and the 
last column to the strongest. 

The QR decomposition at step 50 gives the 
expression: 

S = QT X (13) 

wherein Q is a 2N x 2N unitary matrix, and Ti is a 2N x K 

upper-triangular matrix having the structure T x = I , in 

which T is a K x K upper-triangular matrix. Q can be 
partitioned into Q = [q x Q 2 ], wherein Q x is a 2N x K 
matrix having a column span equal to the column span of 
S, and Q 2 is its complement. In practice, explicit 
computation of Q is preferably avoided, by applying the 
Givens rotations to both S and x . 

Substituting the representation of equation (13) 
into equation (8), we are left with the following model: 

z = Tb + v-l (14) 

wherein z := Q^x , and v x = Q^n . Thus, the transformed 
noise vector Vi remains Gaussian and uncorrelated in 
character . 
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Because T is upper-triangular, the last component of 
b is determined simply, at a last element determination 
step 54, by: 



Assuming that the columns of S were properly ordered at 
step 50, this should give the symbol value for one of the 
strongest user signals, so that the likelihood of error 
10 in determining b k is relatively small. 

The decision value found at step 54 is substituted 
back into equation (14) to obtain a dimensionally-reduced 
model, at a row and column elimination step 56: 



Z (K-1) = T (K-l) b (X-l) + V (K-1) (16) 



Here fcr ' contains the first K-l components of b, T (K " ' 
contains the upper-left (K-l) x (K-l) components of T, and 
the K-l components of z (K_1) are given by 
zf~ 1] = zf ] - T iK b K . The model of equation (16) is again 
triangular in structure. Thus, it is now possible to 
repeat steps 54 and 56, in order to find and eliminate 
b k _ x . The method then proceeds iteratively in this 
manner until all of the elements of b have been found. 

The decision feedback method of Fig. 4 can also be 
applied, mutatis mutandis, to the asynchronous case 
discussed above. For large observation windows (large 
M) , the matrix S is a sparse, banded matrix, having at 
most N+L non-zero entries in each column, located near 
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the matrix diagonal. Banded structures can be exploited 
to reduce the computational cost of QR decomposition. S 
as defined above, however, loses its banded structure. 
The banding can be restored if x is redefined in such a 
5 way that the real and imaginary parts of successive 
samples are interleaved: 

x := ((x^, Ui) 3 , CxaJg,, (x 2 ) 3 , . . . , (x N ) m , (x N )zY (17) 

10 Now the first row of S will contain the real part of the 
first row of S, the second row of S will contain the 
imaginary part of the first row of S, and so forth. The 
real model of equation (8) is still preserved by this new 
structure, so that the methods of decorrelat ion and 

15 decision feedback detection described above can be 
applied here, too. 

Although preferred embodiments are described 
hereinabove with reference to multi-user detection of 
wireless DS-CDMA signals, with BPSK modulation and 

20 complex PN code waveforms, the principles of the present 
invention may be applied to other multi-user detection 
problems. For example, when base station 22 uses 

multiple antennas for spatial diversity, the individual 
user signatures will be complex even when the PN 

25 waveforms are real-valued, due to the inherent complex 
response of the antennas. 

It will thus be appreciated that the preferred 
embodiments described above are cited by way of example, 
and that the present invention is not limited to what has 

30 been particularly shown and described hereinabove. 
Rather, the scope of the present invention includes both 
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combinations and subcombinations of the various features 
described hereinabove, as well as variations and 
modifications thereof which would occur to persons 
skilled in the art upon reading the foregoing description 
and which are not disclosed in the prior art. 
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